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Abstract — Natural draught cooling towers are widely used in modern thermal and nuclear power plants. Because of their 


shear size and susceptibility to horizontal stresses, these towers with very thin shell thickness are unique constructions. The 


boundary conditions should be regarded as free at the top and fixed at the bottom. The cooling tower's material 
characteristics include a young modulus of 31GPa and a Poisson Ratio of 0.15. Finite Element Analysis was used to 


analyse these cooling towers for seismic and wind stresses. Wind loads on these cooling towers have been estimated in the 
form of pressures using the design wind pressure coefficients as supplied in IS: 11504-1985 code together with the design 
wind pressures at different levels as specified in IS: 875 (Part 3) - 1987 code. Ansys 18.2 was used to conduct the analysis. 


The study yields maximum deflection and maximum equivalent stress. 
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I- INTRODUCTION 
The natural draught cooling tower is a critical 
component in both thermal and nuclear power plants. 
These are massive buildings with thin shell structures. 
Cooling towers are susceptible to their own weight as 
well as dynamic loads such as seismic motion and wind 
impacts. In the absence of seismic loading, wind is the 
primary loading for natural draught cooling tower 
design. There has been a lot of study work described in 
the literature on the seismic and wind load on cooling 
towers [1—5]. G. Murali et al., [1] Wind load response of 
a cooling tower. He examined the two cooling towers, 
which were 122m and 200m above ground level, 
respectively. They computed quantities such as 
meridional forces and bending moments. D.Makovika 
investigated the Response Analysis of an RC Cooling 
Tower under Seismic and Windstorm Effects in Acta 
[2]. The 
displacement and internal force envelopes owing to 


Polytechnica estimated values of the 
seismic loading states are compared to the envelopes of 
the loading states due to dead, operational, and live 
loads, wind, and temperature actions. After establishing 
a finite element model, the mechanical properties of the 
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tower under gravity, temperature load, and wind loads 
are examined. A. M. El Ansary [3], A. M. El Ansary [4], 
A. M. El Ans The goal of this work is to create a 
numerical tool capable of attaining an optimum shape 
and design of hyperbolic cooling towers by combining 
an in-house created non-linear finite element model with 
a genetic algorithm optimization approach. The 
influence of asymmetric imperfection on the earth quake 
response of a hyperbolic cooling tower was investigated 
by R.L.Norton [4]. Shailesh S[5], software programme 
used in a real application by taking into account the 
problem of natural draught hyperbolic cooling towers. 
The major goal is to show that the column supports on 
the tower can be replaced by similar shell pieces, 
allowing the software created to be readily used. 


Ii- DESCRIPTION OF THE GEOMETRY OF 
THE TOWER 


The tower's overall height is 200 metres. The base, 
throat, and top radii of the tower are 136 m, 85.27 m, 
and 88.41 m, respectively, with the throat positioned at 
68 percent, 71 percent, 74 percent, 77 percent, and 80 
percent of total height above the base, as shown in Fig. 
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1. It features a 240 mm thick shell wall. Other models' 
size and RCC shell thickness vary in relation to the 
reference tower. The cooling tower's boundary condition 
has been top end free and bottom end fixed.General 
equation of hyperbola is given by 

x? y* 


2 i 


a 
Where x is the horizontal radius at any vertical 
coordinate, y is the origin of coordinates given by the 
centre of the tower throat, a is the radius of the throat, 
and b is any hyperboloid characteristic dimension. 
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Hit- EARTHQUAKE FORCES 


The seismic analysis will be performed in line with IS: 
1893 for lg (g: Gravity acceleration 9810 KN/m2). The 
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seismic study of the shell will be carried out using the 
response spectrum approach and the modal analysis of 
the hyperbolic cooling towers. The response spectrum 
approach will be used for earthquake study of the fill 
supporting structures (RCC frames). The following 
factors were considered for the calculation of the Design 
Spectrum in accordance with IS 1893(Part I)-2002. 


IV- WIND LOADS 


The wind pressure at a particular height [Pz] will be 
calculated in accordance with IS: 875 (part 3)-1987. The 
basic wind speed (Vb) at 9.2m above mean ground level 
will be used to calculate the design wind pressure at a 
specific height. The risk coefficientK1=1.06 will be used 
to compute the design wind speed (Vz) at a height z. For 
coefficient K2, terrain category 2 from IS:875 (part-3)- 
1987 will be used. For design purposes, the wind 
direction will be the one that results in the worst load 
condition. For the tower under consideration, K3 will be 
1. 


V - FINITE ELEMENT MODELING 


Finite element analysis is used due to the complexity of 
the material characteristics, boundary conditions, and 
tower construction. The cooling towers’ finite element 
analysis was performed with ANSYS 18.2. The study 
was carried out with the help of ANSYS workbench. In 
this work, just the shell section of the cooling towers 
were modelled, and fixity was assumed at the base. 


VI - MATERIAL PROPERTIES FOR ANALYSIS 
OF CT 


> Young modulus: 31Gpa 
> Poisson Ratio: 0.15 


VII- TABULATION & RESULTS 


CT 1= Location of throat at a distance of 68% of total 
height of cooling tower measured from base 


CT 2= Location of throat at a distance of 71% of total 
height of cooling tower measured from base 


CT 3= Location of throat at a distance of 74% of total 
height of cooling tower measured from base 


CT 4= Location of throat at a distance of 77% of total 
height of cooling tower measured from base 
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CT 5= Location of throat at a distance of 80% of total 
height of cooling tower measured from base 


Static Analysis 





Fig: 2 Geometry Fig: 3 Boundary condition 





Fig:5 Deflection in CT1 Fig:6 Equivalent Stress in CT 1 


Table 2: Static analysis 


Total deformation | Equivalent Stress 
(mm) (MPa) 


Type of 
cooling 
tower 


5.6714 
57.149 





VUI - CONCLUSION 


The following is the major goal of cooling tower 
analysis work. In the current work, FEA of ten cooling 


towers, CT1 to CT5, was performed to analyse 
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deformation and equivalent stresses. The _ static 
deformation is least for CT 3 for zone 2 and CT 3 for 
zone 4, i.e. when the neck is positioned at 74% of total 
height measured from the bottom.The Equivalent 
stresses in static analysis 1.e. (self-weight) are observed 
to be less for CT 2 for zone 2 and CT 1 for zone 4. In the 
free vibration study, it was discovered that CTS has the 
least deformation when compared to others in zone 2 and 
CTS when compared to others in zone 4. According to 
the seismic study, the equivalent stress is lowest for CT 


2 in zone 2 and highest for CT 2 in zone 4. 
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